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Searches for tetraquarks and mesonic molecules in lattice QCD are 
briefly reviewed. In the light quark sector the most serious candidates are 
the lightest scalar resonances a, k, aq and /q. In the hidden-charm sector 
I discuss lattice simulations of X(3872), y(4260), r(4140) and Z+(4430). 
The most serious challenge in all these lattice studies is the presence of 
scattering states in addition to possible tetraquark/molecular states. The 
available methods for distinguishing both are reviewed and the main con- 
clusions of the simulations are presented. 

1. Introduction 

Some of the observed resonances, i.e. light scalars [1] and some hidden- 
charm resonances [2] , are strong candidates for tetraquarks [qq\ [qq] or mesonic 
molecules {qq){qq)- Current lattice methods do not distinguish between 
both types, so a common name "tetraquarks" will be often used to denote 
both types of qqqq Fock components below. 

In order to extract the information about tetraquark states, lattice QCD 
simulations evaluate correlation functions on x T lattice with tetraquark 
interpolators O ~ qqqq at the source and the sink 

Q,(t) = (0|O.(t)O](0)|0)p^^/-^ J^ZrZf e-^"* n = l,2,... (1) 

n 

If the correlation matrix is calculated for a number of interpolators C'j=i,..,Ar 
with given quantum numbers, the energies of the few lowest physical states 
En and the corresponding couplings Zf = (OlOjln) can be extracted from 
the eigenvalues A"(t) = e-^"(*-*o) g^Yid eigenvectors of the generalized 

eigenvalue problem C(i)iP(f) = A"(t, to)C'(*o)^(*)i as discussed in [3j. 
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In addition to possible tetraquarks, also the two-meson scattering states 
M1M2 unavoidably contribute to the correlation function and this presents 
the main obstacle in extracting the information about tetraquarks. The 
scattering states Mi{k)M2{—k) at total momentum p = have discrete 
energy levels 



in the non-interacting approximation when periodic boundary conditions in 
space are employed. 

The resonance manifests itself on the lattice as a state in addition to the 
discrete tower of scattering states ([2]) [4,^|6j and it is often above the lowest 
scattering state (at E ~ Mi + M2 for S-wave decay). So the extraction of a 
few states in addition to the ground state may be crucial. However, many 
simulations extract only the ground state energy Ei using a conventional 
exponential fit {0\O {t)0^ {0)\0) oc e'^^^ at large t. 

Once the physical states are obtained, one needs to determine whether 
a certain state corresponds to a one-particle (tetraquark) or a two-particle 
(scattering) state and the available methods to distinguish both are reviewed 
in the next Section. 

2. Methods to distinguish one-particle and scattering states 

I am listing the available methods, which may be complementary: 

• For a one-particle state n the coupling is expected to be almost 
independent of the lattice size L, i.e. Z"(Li)/Zf (L2) ~ 1. For a two- 
particle state n = M1M2 one expects (Li)/Zf (L2) ~ (L2/Li)3/2 if 
the range of interaction between Mi and M2 is much smaller than L 
[H El m El [8] . But this method leads to a reliable distinction only in 
presence of long stable plateaus, as cautioned in [9|. 

• One can distinguish whether the ground state is a one-particle or a 
two-particle state from the time-dependence of the Cjj(t) near t ~ 
T/2 at finite temporal extent T. In case of (anti)periodic boundary 
conditions Cii{t) = \Zl\'^ [e~^^^ + {t -^T — t}] for one-particle ground 
state and Cii{t) = \Z}\^ e--^i* + |Z/|2 e-'"Mi*e~"'A/2(T-t) ^ |^ ^ 

{Zf = {M\\Oi\M2)) for two-particle ground state [8j. The criteria for 
distinguishing the excited states is discussed in [5j . 

• An attractive interaction between two particles in a scattering state 
manifests itself by a scattering length a > 0. A formation of a bound 
state below certain m-,^ can be identified by the change of sign of a from 
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positive to negative as m-,^ is lowered [lOl lllj. The scattering length a 
for S-wave scattering M1M2 can be determined from the energy shift 
/S.E = El — niMi — on a finite lattice [11] . 

• Certain non-conventional spatial boundary conditions have specified 
effects on the one- and two-particle energies, which allow to distinguish 
both types [H]. 

3. Light scalar resonances 

It is still not established whether the lightest scalar mesons a, k, ao(980) 
and /o(980) are conventional qq states or have an important qqqq Fock com- 
ponent, as strongly supported by some phenomenological studies [1|. The 
tetraquark interpretation implies that the 1 = 1 state (ussd) is heavier than 
the 7 = 1/2 state (udds) due to rus > md, in agreement with experimental 
ordering ?Ti.ao(980) > "^k- One the other hand, the conventional ud and us 
states can hardly explain the observed mass ordering. The tetraquark in- 
terpretation also naturally explains the large observed coupling of ao(980) 
and /o(980) to KK, which is due to the additional valence pair ss. 

All lattice simulations that look for tetraquark Fock component of light 
scalar mesons are quenched except for [51 1^. All take tetraquark source/sink 
and omit the disconnected contractions in order to look for states with four 
valence quarks. The disconnected diagrams are omitted also since they are 
expensive for numerical evaluation and since they are often noisy. 

• Prelovsek et al. performed the Nf = 2 dynamical and quenched simu- 
lation and extracted three lowest energy states in non-exotic / = 0, 1/2 
and the exotic / = 2,3/2 channels using the variational method and 
a number of and {qq){qq) interpolators [H [5]. The ground 
state in all channels is found to be the scattering state Mi(0)M2(0) 
(M1M2 = TTTT or Ktt), as demonstrated using the time dependence 
of the diagonal correlators. The resulting Zf^L) is also roughly con- 
sistent with the expectation for a scattering state ^^^(12)/^^^ (16) ~ 
(16/12)3/2 (the ratios for 7 = 0, 1/2 have sizable errors, which do not 
allow to make a distinction) [5j. One of the states in all the channels 
is close to Mi(^)M2(-^) state. 

Additional light states are found in / = and 7 = 1/2 channels, which 
may be related to observed a and k resonances with strong tetraquark 
components. The mass dependence of these candidates for a/n on m,r 
are in qualitative agreement with prediction of unitarized ChPT |13j . 
A simulation which takes into account also the disconnected diagrams 
will be needed to verify whether the additional states in 7 = 0,1/2 
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channels are not some kind of unknown artifacts related to the omis- 
sion of disconnected diagrams. In the repulsive / = 2,3/2 channels 
no light state in addition to the scattering states Mi(0)M2(0) and 
Mi(^)M2(— ^) is found, which is consistent with no experimentally 
observed resonances in these two channels. 

• Mathur et al. extract three lowest states in / = channel from a single 
TTTT correlator using the sequential Bayes method [6] and a quenched 
simulation. The ground state energy is consistent with 7r(0)7r(0) and 

its coupling is consistent with scattering state Z^(12)/Z^(16) ~ (16/12)^/^. 
The energy of the third state is consistent with 7r(^)7r(— ^). They 
find an additional state in between, which behaves according to a one- 
particle expectation Z^(12)/Z^(16) ~ 1. This state is a candidate for 
the observed a resonance with a strong tetraquark component. The 
presence of an additional state needs to be confirmed by a simulation 
that takes into account the disconnected contractions. 

• Suganuma et al. extract the ground state from a single [qq] [qq] corre- 
lator in / = channel [13]. They employ the conventional and hybrid 
boundary conditions which indicate that their ground state is a tttt 
scattering state. 

• Alford and Jaffe extract the I = 0, 2 ground state energy Ei{L) from 
a vrvr correlator for a number of lattice sizes L |15j. They argue that 
E[^'^{L) is in accordance with a scattering state, while E{^^{L) de- 
parts from the expected behavior for scattering states and may be an 
indication for a. 

4. Hidden charm resonances 

The most prominent tetraquark candidate is the charged Z"'"(4430) res- 
onance, discovered by Belle [IB]: it decays to ■n~^ij)\ so it must have a min- 
imal quark content ducc, but it has not been confirmed by Babar |17| . 
I will also discuss the observed neutral hidden charmonium resonances 
X(3872), y(4260) and y(4140) [21112], which are candidates for tetraquarks 
or mesonic molecules, although here the charmonium cc Fock component can 
not be straightforwardly excluded based on the charge alone. 

• Chiu and Hsieh [7] simulated states cgcg, cscs, cccc and cqcs {q = n, d) 
with J^^ = V^~^ and J^^ = 1 . They used quenched simula- 
tion with overlap valence quarks and omit the disconnected diagrams. 
They extracted only the ground state energy Ei and coupling Z^{L) 
from diagonal correlator Cii{t) at two different L = 20,24. 
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The ground cqcq state with J^'^ = 1++ is found at 3890 ± 30 MeV, 
which is indeed close to the mass of the X (3872) . They find (20) /Z^ (24) ~ 
1, indicating a one-particle (tetraquark/molecular) state. Note how- 
ever, that the lowest DD* S-wave scattering state with E = mo -\- 
rriD* ~ 3879 MeV is extremely close and that it should be found 
in addition to the one-particle state before the indication for the 
tetraquarks/molecules can be fully trusted. 

The cscs state with J'^*-' = l"*"*" was found (predicted) at 4100 ± 
50 MeV and (20)7^^(24) ~ 1. A state with similar properties 
y(4140) was indeed later observed by CDF [18]. The state is again 
very close to the scattering threshold + nT-j/^ — 4117 MeV, so the 
scattering state has to be found also in order to trust the existence of 
tetraquark / molecule. 

The ground cqcq state with J^^ = 1— is found at 4238 ± 31 MeV, 
which is indeed close to mass of the Y (4260) . They find (20) /Z^ (24) ~ 
1, indicating a one-particle (tetraquark/molecular) state. 

In all three cases a one-particle nature was deduced from (20) /Z^ (24) ~ 
1. However, the cautionary remarks concerning Z(L) ^ have to be 
kept in mind before concluding that tetraquark/molecule really exist. 

• The quenched simulation [11] was to my knowledge the only one aimed 
at the very interesting state Z"'"(4430). The quantum numbers of these 
state are not established experimentally, but since it is very close to the 
DiD* threshold, the simulation [llj is carried out in = 0~, 1~, 2~ 
channels. The scattering lengths a are extracted with the help of 
asymmetric box Li x L2 x L3 x T, which allows for a variety of spatial 
momenta ki = The most reliable results are obtained for = 0~, 
where the attractive interaction between Di and D* is found and 
a > 0. But the a does not change sign with falling m.,^ and the 
authors conclude that the attraction is probably to weak to form a 
loosely bound state. 

• Liuming Liu determined a number of S-wave scattering lengths a for 
scattering between heavy-light, heavy-heavy and light-light mesons 
using 2 + 1 dynamical simulation [10] . She studies only channels where 
no disconnected diagrams are present. The scattering lengths are 
determined from energy shifts AE = Ei — mjvfi — rnM2 ^he ground 
state energies Ei are obtained using the M1M2 interpolators. As far 
as tetraquarks/molecules are considered, the most interesting result 
comes from the D'^D^* channel with 1=1, where a changes sign 
at m^r — 280 MeV. This may be an indication for an existence of a 
loosely bound state at < 280 MeV. 
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• The 2 + 1 dynamical simulation of Ehmann and Bali [19j was actually 
not aimed at searching for tetraquarks/molecules but to study the 
mixing between charmonia {J ijc, Xc) and DD (D stands for D, D*, 
Di) states. They compute the a full correlation matrix with cTc as well 
as {cTq){qrc) interpolators in J^'-" = '",1 ,2"'"+ channels, taking 
into account also all disconnected contractions. Using the variational 
method they determine the spectrum and also the components cTc 
and {cTq){qTc) of the charmonia and DD physical eigenstates. They 
do not find any states in addition to the expected charmonia and 
scattering states and they do not attempt to establish whether their 
resulting states are one-particle or scattering states. Let me note that 
a nonzero coupling {0\{cTq){qTc)\J/^) , for example, does not mean 
that J/ip has a sizable tetraquark component, since {cTq){qTc) and 
cTc Fock components mix via singly disconnected contractions in |19j . 

5. Conclusions 

Proving a sizable tetraquark or molecular Fock component in a hadronic 
resonance using lattice QCD simulation is not an easy task. A resonance 
appears as a state in addition to the discrete tower of scattering states. So 
the extraction of few states in addition to the ground state is expected to 
be crucial. Given the resulting physical eigenstates, one needs to determine 
whether a certain state corresponds to a one-particle (tetraquark/molecular) 
or a two-particle (scattering) state, and the available methods to distinguish 
both are reviewed. 

There are some indications for an additional state in / = 0, 1/2 light 
scalar channels, which might correspond to observed a and k with strong 
tetraquark components ^ 13 |^. The corresponding simulations omitted 
the disconnected contractions when calculating correlators with tetraquark 
interpolators in order to study genuine tetraquark states with four valence 
quarks. It would be valuable to verify in the future whether the additional 
states are present also in a simulation which takes into account the discon- 
nected contractions. 

There have been surprisingly few lattice simulation of very interesting 
exotic XYZ resonances, discovered recently in B-factories. Most of simula- 
tions extract only the ground state in a given channel and then try to deter- 
mine whether it corresponds to a one-particle (tetraquark/molecular) state 
or to a scattering state. There is some indication that X(3872), y(4260) and 
y(4140) are tetraquark/molecular states, but future simulation are needed 
to verify that. 
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